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ABSTRACT 

We present observations of the HCN and HCO+ J=l-0 transitions in the center of the nearby spiral 
galaxy NGC 6946 made with the BIMA and CARMA interferometers. Using the BIMA SONG CO 
map, we investigate the change in the Incn/Ico and /hco+/^CO integrated intensity ratios as a 
function of radius in the central kiloparsec of the galaxy, and find that they are strongly concentrated 
at the center. We use the 2MASS Ks band image to find the stellar surface density, and then construct 
a map of the hydrostatic midplane pressure. We apply a PDR model to the observed /hcn/-^hco+ 
integrated intensity ratio to calculate the number density of molecular hydrogen in the dense gas 
tracer emitting region, and find that it is roughly constant at 10^ cm~^ across our map. We explore 
two hypotheses for the distribution of the dense gas. If the HCN and HCO"*" emission comes from 
self-gravitating density peaks inside of a less dense gas distribution, there is a linear proportionality 
between the internal velocity dispersion of the dense gas and the size of the density peak. Alternatively, 
the HCN and HCO+ emission could come from dense gas homogeneously distributed throughout the 
center and bound by ambient pressure, similar to what is observed toward the center of the Milky Way. 
We find both of these hypotheses to be plausible. We fit the relationships between /hcn, ^hco+, and 
/co • Correlations between the hydrostatic midplane pressure and /hcn and /hco + are dem onstrated, 
and p ower law fits are provided. We confirm the validity of a relation found bv lBlitz fc R osolowskvl 
(|2006f ) between pressure and the molecular to atomic gas ratio in the high hydrostatic midplane 
pressure regime (10^-10* cm^"^ K). 

Subject headings: Instrumentation: interferometers — ISM: molecules — galaxies: ISM — galaxies: 
nuclei — radio lines: galaxies 
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1. INTRODUCTION 

How is molecular gas organized in galaxies? In the 
disk of the Milky Way, molecular gas as traced by CO 
is concentrated in giant mole cular cloud s (GMCs) with 
a characteristic size of 50 pc (|Blitzjll987D . These clouds 
have masses aro und IQ-^-lO^ Mq and contain 80% of the 
molecular mass ()Blitj [l993l ) . The hydrostatic pressure 
from the self-gravity of a GMC this size is ~ 10^ cm~'^ K 
{Blitz 1993), which is an order of m agnitu de larger than 
the mean pressure of the local ISM ([Bloem en 1987); this 
leads to the conclusion that GMCs in the Milky Way disk 
are self-Rra .vitating. GMCs in M31 (Rosolowsky 2007) 
and M33 ()Rosolowskv fc Blitd l2004l : IRo s olowsk v et all 
l2007f) have similar properties to those observed in the 
Milky Way 

The densest parts of GMCs are traced by molecules 
with larger dipole moments than CO, such as HCN, 
HCO+, and CS. These dense gas tracers need a molec- 
ular hydrogen density near 10^ cm~^ to be excited col- 
lisionally, 100 times greater than CO. In the disk, emis- 
sion from these molecules is limited to the star-forming 
cores of GMCs (iPlmneiEaD [IMl iWi fc Evan j l2fl03f 
iHelfer fc Blitzlll997an . 

However, the molecular gas composition of the central 
kiloparsec of a spiral galaxy is very different from that 
of its disk. In the Milky Way, widespread CO emission 
in the inner kpc makes it difficult to differentiate be- 
tween a distribution of discrete clouds or one in which the 
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mole cular gas is organized homogeneously (jPame et al.l 
[19871) . The emission from HCN and CS within ~250 
pc o f the Galactic center is distributed sim ilarly to the 
CO (|,Tackson et al.l[l996l : iBaflv et aLlflQSl . suggesting 
that high density regions are common, at least com- 
pared to the average ISM density of ~1 cm~^ in the 
solar neighborhood. In other galaxies, dense gas tracers 
are more concentrated than CO in the central kpc and 
have lower intensities elsewhere ()Gao fc SolomonlT2004al : 
iHelfer fc Blitd 11997 b') , consistent with observations of 
the Milky Way (jHelfer fc Bhtz 1997a). Furthermore, 
the centers of spiral galaxies tend to have high den- 
sities and pressures; bulges have hydrostatic midplane 
pressures 2- 3 orders of mag nitude higher than those of 
their disks (jSpergel fc Blitzlll992f) . iBlitz fc Rosolowskvl 
{2006) (hereafter BROS) found a simple relationship in 
disk galaxies between the hydrostatic midplane pressure 
and the ratio of molecular to atomic gas that holds in 
disks and bulges. 

In this paper, we use new HCN and HCO+ ob- 
servations of the central region of NGC 6946 to 
study molecular gas organization in this nearby spiral 
galaxy. Optically, thi s galaxy is classified as SAB(rs)cd 
(Ide Vaucouleurslll963l), it ha s a moderate starburst in its 
nucleus ("Turne r fc Holll983h . and it is at a distance of 5.5 
Mpc (Tully 198i). At this distance, 1" corresponds to 27 
pc. Observations of NGC 6946 in the R, I, and K bands 
reveal a sma ll bulge (radius 15") within a larger bar 
(radius 63") ()Das et al.ll2003f ). Recent studies revealed 
the ric h details of its nuclear structure; ISchinnerer et al.l 
j200(l ) used high resolution images of CO to map the 
molecular gas spiral in the inner 10" , likely caused by an 
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inner stellar bar 15" i n le ngth (see also NIR i mages in 
lElmegreen et al.lll998l and lKnapen et al.ll2003f ). A fur- 
ther study used HCN emission in the central 2" as a 
proxy for star formation in young clusters still embedded 
in their parent clouds, and combined these data with ob- 
servations of HII regions to map th e total amount of star 
formation (jSchinnerer et al.1l2007[ ). They concluded the 
nuclear star formation is being fed by an accumulation 
of dense gas driven by the inner stellar bar. 

Using HCN and HCO+ observations in conjunction 
with CO and near-infrared data, we investigate the distri- 
bution of dense gas in the central kiloparsec of NGC 6946. 
After reducing our observations, we calculate HCN/CO 
and HCO+/CO integrated intensity ratios, and we es- 
timate the pressure using two different methods. We 
then describe two hypotheses for the distribution of the 
dense molecular gas, and explore the relationship be- 
tween dense gas tracers and hydrostatic midplane pres- 
sure. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. New Observations 

2.1.1. BIMA HCN and HCO+ Observations 

We observed the HCN J = 1 — transition at 
V = 88.631847 GHz simultaneously with the HCO+ 
J = 1-0 transition &t v = 89.188518 GHz us- 
ing the 10-element Berkeley- Illinois-Mary land Associa- 
tion (BIMA) interferometer (| Welch et al.lll996t ) at Hat 
Creek, CA. The observations consisted of five short 
tracks during September and October 2002, totaling 
20 hours of observations in the C and D configura- 
tions. Baselines ranged from 1.9 to 26 kA. The 
pointing and phase center of the observations was 
a(J2000)=20''34™52!33, (5(J2000)=-|-60°09'14'.'2, and we 
used a seven-point hexagonal mosaic with a Nyquist 
spacing of 57" . The correlator was centered near 52 km 
s~^, the LSR velocity of NGC 6946, and configured to 
achieve a spectral resolution of 5.3 km s~^ over a velocity 
range of -120 < w < 210 km s^^ for both the HCN and 
HCO+ observations, which extends beyond the detected 
velo city range of the C O emission (^ —50 <v< 180 km 
s-^ HHelfer etm[2003l) . 

2.1.2. CARMA HCN Observations 

We observed the HCN J = 1 — transition using the 
Combined Array for Research in Millimeter-wave Astron- 
omy (CARMA)^ at Cedar Flat, CA, for 11 hours over 
three tracks in December 2006 and March 2007. CARMA 
is the merger of nine 6.1 m BIMA telescopes with six 10.4 
m Owens Valley Radio Observatory dishes. In Decem- 
ber, the array was in the first light configuration, which 
included baselines as short as 3.3 kA and as long as 31 
kA. In March, the telescopes were in C array, with base- 
lines ranging from 8 kA to 101 kA. We configured the 
correlator with two overlapping 62 MHz windows and a 
resolution of 0.977 MHz (~3 km s~^). This setup corre- 
sponded to a velocity coverage of —140 < w < 230 km 

^ Support for CARMA construction was derived from the states 
of California, Illinois, and Maryland, the Gordon and Betty Moore 
Foundation, the Eileen and Kenneth Norris Foundation, the Cal- 
tech Associates, and the National Science Foundation. Ongoing 
CARMA development and operations are supported by the Na- 
tional Science Foundation under a cooperative agreement, and by 
the CARMA partner universities. 
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Fig. 1. — HCN channel map from combined BIMA and CARMA 
observations. The velocity is marked in the upper right corner of 
each plot in km s~^. Contour levels are drawn at 2,4,6,8. . .times 
the rms channel noise. A filled grey contour is marked at -2 times 
the rms channel noise. The beam is shown in the lower left hand 
corner of the -50 km channel. 

s~^. We used a seven-point hexagonal mosaic with a 30" 
Nyquist spacing of the 10.4 m primary beams. We ob- 
served Uranus and 3C279 to calibrate the flux scale and 
the frequency response of the passband. Approximately 
every 25 minutes throughout each track, we observed a 
nearby quasar to calibrate the time dependence of the 
phase. 

2.1.3. Inter ferometric Data Reduction 

We reduced the H CN and HCO+ data using the 
MIRIAD package (Sa ult et all Il995f ). We combined 
the calibrated M tf- visi bilities and then made robustly 
weighted (Briggi ri995[ ) channel maps with velocity bins 
of width 10 km s~^. We deconv olved the image cube 
using an SDI CLEAN algorithm (|Steer et al.lll984[ ): re- 
stored channel maps from this procedure are shown in 
Figs. [Hand m The synthesized beam for the HCN map 
was 4.7" by 3.7" with a position angle of 61°; for the 
HCO+ map the beam was 8.6" by 7.5" with a position 
angle of 87° . We calculated the noise level by taking the 
average of the rms variation measured in the central ar- 
cminute box over several line-free channels on either side 
of the HCN or HCO+ emission. The rms per 10 km s~^ 
channel was 0.012 Jy/beam for the HCN map and 0.020 
Jy/beam for HCO+ map. Based on our flux calibration 
and comparison of the fluxes from individual tracks, we 
assigned a la error in the flux calibration of ±20%. 

2.2. CO Observations 

Observations of the J = 1 — line of CO dX v = 
115.2712 GHz in NGC 6946 were taken from the BIMA 
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Fig. 2. — HCO"*" channel map from BIMA observations. The 
velocity is marked in the upper right corner of each plot in km 
s~^. Contour levels are drawn at 2,4,6,8. . .times the rms channel 
noise. A filled grey contour is marked at -2 times the rms channel 
noise. The beam is shown in the lower right hand corner of the -50 
km s~^ channel. 

Survey of Nearby Galaxies (SONG) (|Regaii et al.l[200ll : 
iHelfer et al.ll2003l ). These observations have a resolution 
of 1.56 MHz (~4 km s^"'^) over a total bandwidth of 368 
MHz (960 km s^^). The interferometer data was com- 
bined with observations from the NRAO 12m telescope 
at Kitt Peak, h!L? The synthesized beam was 6.0" by 
4.9" with a position angle of 14°. The \a flux calibra- 
tion error was ±15%. 

2.3. Integrated Intensity Maps 

To construct integrated intensity maps of the dense 
gas tracers while adding in as little noise as possible, we 
masked out the emission in regions where the signal is 
below the noise threshold. We started with the channel 
maps Fx{a, 6, v) where X denotes the species and v the 
velocity channel. We constructed independent masks for 
HCN and HCO+ by first creating blurred maps through 
convolution with a 20" Gaussian, and then averaging the 
rms noise in several channels off the line to calculate 
fx.biur- For each of these two species we defined a mask: 

Cx{a, S,v)^l where Fx > 3crx,biur 

= elsewhere. (1) 

We then made integrated intensity maps: 

/x(a,<5) (2) 

V 

^ The National Radio Astronomy Observatory is operated by the 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 



30 20 10 -10 -20 -30 
RA offset (arcsec) 
Fig. 3. — HCN integrated intensity map. Rod/orange colors 
are CO, black contour lines are HCN. Contours are logarithmi- 
cally spaced in both CO and HCN, such that a contour repre- 
sents 2,2\/2,4,4\/2,8. . . times 0.636 Jy/beam km s"! for HCN and 
2,4,8,16. . .times 3.01 Jy/beam km s~^ for CO. These factors are 
the Qx noise levels for each species for a point where all the chan- 
nels on the line are included (see i|3.2l l. The HCN beam is shown 
in the upper left hand corner, and the CO beam is in the upper 
right. 

where dv is the channel velocity width. This masking 
technique preserves noise statistics, but introduces a bias 
in favor of low-level diffuse emission distributed simi- 
larly to the bright emission and against low-level compact 
emission distributed differen tly t han the brightest e mis- 
sion (see lRegan et al.l (|200lD and 'He ifer et all (|200 3^ for 
more thorough discussions). Fig. [3] shows /hcn contours 
plotted over Iqo color, and Fig. 2] plots /hco+ contours 
over Ico color. For these figures, no mask was applied 
to the CO data. 

2.4. 2MASS Image Processing 

We obtained Ks (2.2 ^m) band images of NGC 6946 
from the 2MASS Large Galaxy Atlas (jJarrett et al.l 
l2003l ). Galaxies in this survey have a typical angular 
resolution of 3", with a 1" pixel grid. First, we applied a 
median filter to the image to remove foreground stellar 
contamination. Next, we corrected the im a.ge for dust ex- 
tinct ion from the Milky Way (^^ =0.125) (jSchlegel et all 
mH) as given in NED^. Dust in NGC 6946 will also af- 
fect the observed Ks band image, but dust ext inction 
at the se wavelengths is generally small. Indeed, iReganl 
(|2000f ) calculated a dust optical depth map for NGC 6946 
by using optical and near-infrared images as inputs to 

* The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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Fig. 4. — HCO+ integrated intensity map. Red/orange colors 
are CO, black contour lines are HCO"*". Contours are logarith- 
mically spaced in both CO and HCO"*", such that a contour rep- 
resents 2,2^2,4,4^2,8. .. times 0.982 Jy/beam km s"! for HCO+ 
and 2,4,8,16. . .times 3.01 Jy/beam km s"! for CO. The HCO+ 
beam is shown in the upper left hand corner, and the CO beam is 
in the upper right. 

a radiative transfer model. Though his model did not 
converge in the inner few arcseconds, he found V band 
extinctions peaking near 5 magnitudes just outside that 
region. The dust extinction varied strongly as a function 
of distance from the galactic center. Using the relation 
between the ma gnitude of dust absor ption in V and K for 
the Milky Way (|Schlegel et al.|[l998f ). this corresponds to 
Ak ^0-5 which is about a factor of 2 in flux. 

3. ANALYSIS 
3.1. Comparison to Previous Observations 

Figs. [T] and [2] show HCN and HCO+ diffuse emission 
starting from —20 km s~^ in the north and shifting south- 
ward with higher velocity before falling off near 120 km 
s~^. CO channel maps show similar kinematic structure, 
albeit over a la r ger s patial region (Hclfcr ct al. 2003; 
ISchinnerer eTall l200l . The HCN emission is largely 
confined to the inner most 20" , consistent w ith previous 
maps of this galaxy (|Helfer fc Blitd Il997bl) . Observa- 
tions with a beamize of ~1" show that the HCN emis- 
sion breaks down into a central clump elongated in the 
east-west direction, and a wea ker emitting region 3 " east 
and 1" south of the center ^Schinnerer et al.]l2007t l. 

The integrated intensities we measured are br oadly 
consistent with those compiled bv lMeiierinkI (120061 ) con- 
volved to a single dish 21" beam pointed at the center of 
NGC 6946; he found /rcn = 11 K km s^^, /hco+ = 13 
K km s^^, and Iqo — 220 K km s^^. At this resolu- 
tion, we measured /hcn = 13 K km s^^, /hco+ = 14 K 
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Fig. 5. — Emission intensity ratio /hcn/^CO on the sky coordi- 
nates. The contour levels are shown in the colorbar on the right, 
and the beam is shown in the upper left hand corner. For ease of 
reading, contours starting at 0.04 are also marked with a solid line. 
The features seen at locations away from the central emission that 
are not seen in the HCN map in Fig. |3] have a low value of /hcn. 



km 



and Ico= 280 K km s ^. Differences between 



these two sets of observations can be explained by er- 
rors in flux calibration and single dish pointing. When 
measured with a single-dish telescope, a typical spiral 
alaxy has an IncN/Ico ratio of 0.02-0.06 in its bulge 
Heifer fc Blitlfiggl ICaofc Solomonll2004aD . The inte- 
grated intensities for NGC 6946 result in ratios that are 
well within this range. 

3.2. Ratio Maps 

The Ihcn/Ico ratio is useful because it is a qualita- 
tive indicator of the density of the emitting gas. In order 
to construct ratios of the HCN and CO maps, we con- 
volved the HCN channel map to a 6.0" by 4.9" Gaussian 
beam to match the CO observations, and then created a 
mask based on the CO data as described in ^2.'6l call this 
mask Ceo- We recalculated the HCN, HCO+, and CO 
integrated intensity maps using Ceo in eqn. [^J for the 
remainder of the paper, Ix will refer to the integrated 
intensity calculate with this new Ceo mask. We applied 
the same mask to all three molecules to avoid biasing the 
ratios. The ratio Ihcn/Igo is shown in Fig. [H We also 
constructed independent noise maps, Qx, for /cOj -^hcNj 
and /hco+ using the following formulation: 



(a, 5) 



dv 



/^Cco +iuxlxf (3) 



where dv is the channel velocity width and ux is the 
uncertainty due to flux calibration for each species. 

The ratio in Fig. [5] is plotted only at points where 
^HCN > 2Qhcn and Iqo > 5Qco- The ratio in the 
central 20" peaks at 0.081 ± 0.021 about 1" east of the 
pointing center and has a relatively smooth falloff with 
increasing angular distance from the center; also, the 
three highest contours appear to be elongated along the 
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east-west axis more than along the no rth-south one, as 
expec ted from the high resolution map (jSchinnerer et alJ 
l2007t) . Extensions in the HCN distribution 25" north 
and south of the pointing center are prominent, despite 
appearing weakly in the channel maps. These exten- 
sions have small pockets of /hcn / ^co ratios that peak at 
~0.08, and they coinci de with the rnolecu lar arms seen in 
previous CO maps fi.e. 'Reg an et al.ll2001[ ). Small patches 
of HCN emission further from the center of the map have 
high ratios because of the low intensity of CO emission in 
those regions; they also have relatively little HCN emis- 
sion^ 

iHelfer k, Blitd (|1997bl ) measured the /hcn/^co ratio 
in NGC 6946 and found a central ratio of ~ 0.10 using 
data with nearly identical spatial resolution but some- 
what poorer sensitivity. The structure of their ratio map 
is similar to ours, with the highest ratios along an east- 
west extended structure, but in their observations the 
structure is saddle shaped. It has a peak ratio of 0.19 
and maxima ~5" east and west of the center. These sad- 
dle peaks are not reproduced in our observations, but 
they were the most uncertain features of the earlier map. 

Using the same method, we also created a map of the 
-^HCO+Z-^CO ratio (Fig. [S]). To make this map, we con- 
volved the CO channel map to the larger beamsize of 
the HCO+ observations. The peak ratio in this map is 
0.069 ± 0.018 measured less than an arcsecond from the 
pointing center. The higher peak ratio in the /hcn/^co 
map is due to its smaller beamsize and is not physically 
significant. The extensions to the north and south of the 
pointing center that appeared in the /hcn / ^co raap are 
present here as well. 

To compare the distribution of /hcn and /hco+j we 
constructed a final ratio map of the emission from the 
two dense gas tracers by convolving the HCN map to 
the HCO+ beamsize (Fig. [71). The ratio is only plotted 
at points where /hcn > SQhcn and /hco+ > 2Qhco+- 
Overall, the /hcn//hco+ ratio is close to 1 over the cen- 
tral 20"; since the two molecules have similar critical 
densities this was expected. In the cental regions there 
is a weak trend toward higher ratios in the southeast and 
lower ones in the northwest. In the extensions, the ratio 
appears to fall as low as 0.4, but the relatively low in- 
tensity of the emission from these features leads to large 
uncertainty in the ratio. 

3.3. /hcn//cO; ^hco+AcOj o^nd /hcn//hco+ Radial 

Falloffs 

Although the /hcn //co ratio map is roughly a map of 
the density of the emitting gas, nearby points in the map 
are not independent; their degree of correlation is quan- 
tified by the synthesized beam. To answer the question 
of how this qualitative measure of the density varies with 
radius, we constructed a map with uncorrelated data 
points. We began by bilinearly interpolating the values 
of /hcn , /co , Qhcn and Qco on a hexagonal grid with 
6.0" x4.9" spacing. This resampling results in intensities 
and errors that are only weakly correlated, since they 
are separated by a full beamwidth in the CO map, which 
has the larger beamsize. Throughout this paper, we will 
treat the points on this grid as uncorrelated. 

In order to find how the ratio changes with distance 
from the galactic ce nter, we corr e cted for the inclina- 
tion of the galaxy, i. iRegan et al.l (|2001f ) found i — 54° 
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Fig. 6. — Same as Fig.[5]but for /hco+/^CO- 




-sol \ \ \ \ \ \ \ I I I 

30 20 10 -10 -20 -30 
RA offset (orcsec) 

Fig. 7.— Same as Fig.[5]but for /hcn//hco+- The 0.9, 1.0, and 
1.1 contours are marked with sohd hnes. 

and a position angle of 65°; iTullvl ()1988f ) found that 
the distance to NGC 6946 is 5.5 Mpc. Using these 
parameters, we deprojected the coordinates on the sky 
to physical coordinates centered on the pointing cen- 
ter of the observations. The inclination and position 
angle of NGC 6946 are a subj ect of disagreement, as 
IConsidere fc Athanassoulal 11988) found that NGC 6946 
has i = 34 ± 8° and a position angle of 69±15°. We 
deprojected our data with both sets of parameters, and 
found that our results do not strongly depend on the 
assumed orientation of the galaxy. 

Fig. [H ( Top) plots the /hcn / /co ratio as a function of 
average physical distance from the pointing center for the 
innermost points in the independently spaced hexagonal 
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Fig. 8. — Top: Emission f5iY9isily ratio /hcn/^CO on the 
6.0" x4.9" hexagonal grid. The radial coordinate is the deprojected 
physical radius averaged over a circular beam centered at each grid- 
point. The units of the ratio are K km s~^/K km s~^. The line is 
the fit from eqn.|4] Middle: The same plot for /hco+/^CO on the 
8.6" X 7.5" grid. The line is the fit from eqn. [5] Bottom: The same 
plot for /hcn/^hco+- The line is the fit from eqn. |6] 

grid. The average physical distance, i?avg is determined 
for each point by averaging the radius over an eUipse 
with major and minor axes equal to the grid spacing; 
we did this primarily to get an accurate average distance 
for the central point. For the innermost points the error 
is dominated by the uncertainty in the flux calibration, 
while at larger radii the channel noise is dominant. The 
ratio falls off strongly as a function of i?avg- Ratios at 
radii larger than 380 pc suffer from systematic errors, 
and are not included in this analysis. If we sum the 
integrated intensities of all the points with i?avg between 
380 and 500 pc, we find a ratio of 0.018 ± 0.006. This 
ratio is comparable to that found at the outermost points 
in Fig. [H] ( Top) , which suggests that the ratio flattens out 
at larger i?avg- However, non-axisymmetric features with 
elevated ratios such as the molecular arms are a more 
likely reason for this trend. 

We performed a linear fit on the /hcn / -^co ratio as 
a function of i?avg for the inner 380 pc, and found: 



-^HCN 

Ico 



0.091 ± 0.016 - (2.0 ± 0.6) x 10"'' i?a 



(4) 



The slope of this ratio is different from zero with 3cr sig- 
nificance. Note that by reducing the ratio to a function 
of i?avg , we ignored any differences in physical conditions 
between points in the bar and those outside of it. We 
used a linear model because of its simplicity, but some 



other function lik ely fits the dat a better. 

As noted in iHelfer fc Blitg (|1997b[ ). the central 
^HCN /Ico ratio in NGC 6946 is comparable to that mea- 
sured in the inner 63 pc of the Milky Way, 0.081 ± 0.004 
( Jack son et al.l Il996l ) . This suggests that the physical 
conditions, and therefore the molecular gas distribution, 
in the two systems are similar. 

To make the same plot for the Ihco+Z IcOi we con- 
structed a new hexagonal grid because the HCO"'" data is 
of a lower resolution than the other species. We followed 
the prescription described above and create a new grid 
with 8.6" X 7.5" spacing. The falloff in the /hco-i- /^CO ra- 
tio with radius on this grid is shown in Fig. [8] (Middle). 
The fit to this ratio is: 



iHC04 

^co 



= 0.081 ± 0.015 - (1.1 ± 0.5) X 10"'' R 



avg- 



(5) 



The slope of this ratio is different from zero with 2a 
significance. 

The strong falloff in the /hcn / Ico and /hco+ /Ico ra- 
tios suggest conditions in the inner 100 pc are quite dif- 
ferent from those further from the center, and markedly 
different from those in the disk. Thus, there is no rea- 
son to expect that the molecular gas in the nucleus is 
distributed similarly to that in the disk. 

Finally, on the coarse grid, we plot the ratio 
Ihcn / Ihc'0+ {Fig^ {Bottom)). The fit to these points 
is given by: 



I 



HCN 



0.80 ± 0.26 + (6.9 ± 10.0) x 10"* R 



, , , (6) 

-'HCO+ 

Note the large uncertainty in the slope; it is consistent 
with zero. We performed this last fit mainly as a consis- 
tency check; since HCN and HCO"*" are both dense gas 
tracers with similar critical densities, their ratio should 
be roughly constant with i?avg- This is indeed what we 
observe. 

3.4. Molecular Line Ratios 

Since emission from CO, HCN, and HCO+ all trace 
the molecular gas, there should be corr elations between 
their in tegrated intensities. Following iGao fc SolomonI 
()2004bf l. we fit the relationships between the integrated 
intensities of the three dense gas tracers. We cannot use 
the linear fitting routine from the previous section 
because the points have errors in all parameters; fur- 
thermore, power laws are straightforward to fit in log-log 
space, but there the errors are asymmetric. For these 
reasons, we fit the relationships between the integrated 
intensities using a Monte Carlo fitting routine. A posi- 
tion for each point was randomly chosen using the cal- 
culated errors, and then a linear least squares fit was 
performed to solve for the fit parameters. For each fit, 
we repeated this process lO'^ times and took the mean 
and standard deviation of the fitted parameters. 

We begin with the relationship between /hcn and Ico 
(Fig. El (Top)): 



logio /hcn 



(1.7±0.3)logio/co -3.1±0.8. (7) 

Points in the lower integrated intensity regions have com- 
parably larger errors in /hcni so the poor fit in that 
region is of little significance. The overall fit suggests a 
power law relation with a slope significantly steeper than 
linear. 
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Fig. 9. — Top: /hcn plofr(le(2)aI a function Jf Ico- Error bars for 
both molecules include random noise and the error in flux calibra- 
tion. We plot the fit from eqn. [7] as a solid line. Bottom: /jjCO-l- 
plotted as a function of Iqq . Error bars for both molecules include 
random noise and the error in fiux calibration. We plot the fit from 
eqn. |8]as a solid line. 



The same plot for HCO+ is shown in Fig. [9] {Bottom). 
The fit for this relation is given by: 



login -^HC04 



(0.8±0.2)logio/co -0.9±0.6. (8) 



The single power law appears to fit less well here than in 
Fig. [9] ( Top), but there are not enough independent data 
points to try a more complicated function. 

These fits are an alternate method of presenting the 
variation in the /hcn/^co and /hco+/^co ratios. In 
other words, the value of the slope in eqn. [7] was not a 
surprise; it was clear from the shape of the ratio map cen- 
ter (Fig. [5]) that the slope would be steeper than linear. 
We expected a similar slope in the HCO"'" fit, but the best 
fit to the data was a near-linear power law. Since this 
fit is on the coarser grid, there are fewer points to con- 
strain the relation, and the effect of fluctuations in the 
integrated intensities is greater. We suspect this results 
in an incorrect slope. The relatively poor quality of the 
fit in Fig. [5] (Bottom) lends credence to this argument. 

Physically, what does this analysis of the molecular line 
integrated intensities and ratios imply about the distri- 
bution of molecular gas? A steeper than linear correla- 
tion between a dense gas tracer and CO suggests that the 
proportion of dense gas rises as /co increases, and also 
as distance from the center decreases. To understand 
what could cause this, we now investigate the physical 
conditions in the emitting region. 



3.5. Hydrostatic Midplane Pressure 

We constructed a measure of the hydrostatic midplane 
pressure to compare with the dense gas tracer emission. 
We first convolved the corrected Ks image to the same 
beam as the CO observations, and sampled the resulting 
image on the hexagonal grid. The stellar surface density 
at each point is then given by: 

^°gio * -2 = -^■^^^Ks + 9.62 + logio cos i (9) 
Mq pc 

where fixs is the magnitude of flux in the Ks band, in 
magnitudes per square arcsecond (BR06). 

BR06 used the following formula for the hydrostatic 
pressure, Phyd, in the midplane of a disk: 

0.5 



'^^^=272 cm-3 K 



M, 



PC" 



M, 



o pc" 



km s 



-0.5 



(10) 



The calculation of the gas and stellar surface densities, 
Sg and respectively, is described in the same paper. 
Eqn . II 01 provides a rough estimate of the hydrostatic mid- 
plane pressure, since it was derived for a disk geometry 
and we used it on an elliptical stellar distribution. We 
used the CO integrated intensity to determine the molec- 
ular gas content, despite the possibility that using a con- 
stant X-factor with CO will ov erstate the molecular gas 
conte nt in high density regions (|Reganll200(3l : lRegan et all 
l2001h . We did this because we need a measure of the 
pressure that is independent of /hcn and /hco+- The 
molecular gas surface density is larger than 50 Mq pc~^ 
throughout the inner 50", and within 10" of the pointing 
center it is larger than 500 Mq pc~^ We used a con- 
stant central HI gas surface density (corrected for incli- 
nation) of 3.9 Mq pc"2 (Tacconi & Young 1986), which 
is much smaller than the molecular gas durface density. 
The gas velocity dispersion, Vg is seen in f ace-on galaxies 
to be roughly constant (Shostak & van d er Kruitl[l98l 
iDickev et ailllQQOl : iPetric fc Rupem i200lD : here we as- 
sume Vg is 8 km s~-^ throughout NGC 6946. Note that 
this is generally interpreted as a large-scale turbulent ve- 
locity width, which is different from the velocity width 
within an indi vidual cloud. This dist inction will be im- 
portant below. iRegan fc Vogell (|1995l ) deconvolved a K'- 
band image of the center of NGC 6946 into components, 
and found evidence for a 3.5" oval distortion inside a 15" 
bulge. We used the size of this smaller component to 
estimate the stellar scale height, /i*. 

The hydrostatic midplane pressure is the pressure re- 
quired to balance the gravitational force of the stars; this 
pressure generally comes from turbulent or thermal mo- 
tion of the gas. Note that Phyd is the minimum required 
pressure, and the actual pressure in the midplane could 
be higher. Fig. 1101 shows the hydrostatic midplane pres- 
sure calculated as a function of position on the sky; it 
peaks at ~ 5 X 10^ cm^"^ K and falls off by more than 
an order of magnitude in the inner 10" . The uncertainty 
in the Ks band flux due to dust discussed in t j2.4l cor- 
responds to a small uncertainty in the pressure because 
the stellar surface density enters with a square root. 

3.6. Number Densities from Photon Dominated Region 

Models 
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10 -10 
RA offset (orcsec) 
Fig. 10. — A plot of logj^Q(Pijy(j/fc) as a function of position on 
the 6"x4.9" hexagonal grid. The units of P^yd/fc are cm~^ K. The 
grid points are overlaid in light grey x's. 

We used a code that relies on the /hcn / Iwo+ line ra- 
tio as sampled on the 8.6" x7.5" grid to calculate a num- 
ber density of molecular hydrogen. This number density, 
niine, is specifically indicative of the physical conditions 
in the dense gas tracer emitting region, since it depends 
only on emission from tracers with high critical densities. 

We utilized a Photon Dominated Region (PDR) code 
to calculate the thermal balance (with line transfer) self- 
consistently with the chern ical balance through iteration 
(|Meiierink &: Spaansll2005[) . Both the density and ultra- 
violet radiati on field Gp w ere varied over a large grid 
of models (M eijerink et al.l i2007). Using the chemical 
and thermal structure of the PDR models along with a 
radiation transf e r cod e ([Poelman fc Spaanil2005l [2006D . 
iMeijerink et al.l (|2007t l calculated the line intensities and 
ratios for a number of species including HCN and HCO"'". 
We compared these ratios to those observed for the cen- 
ter of NGC 6946 to find nune- 

We found points in the inner 380 pc have number den- 
sities 5x 10"* cm^'^ < riiino ^ 10^ cm~'^ and average cloud 
temperatures T 20 K for Gq = 100. Choosing larger 
values of Go will raise T; for example, if Go = 1000, 
T ~ 45 K. The v a lue of Go does not affect nuno- Note 
that IStacev et al.l (|1991h constrained the radiation field 
for the central 55" of NGC 6946, and found Go - 500, 
consistent with our assumptions. We plot niinc as a func- 
tion of Phyd in Fig. [11] The number density of the dense 
gas does not appear to vary as a function of Phyd or R„ 
We explore the implications of this below. 




4. DISCUSSION 

We discuss two plausible hypotheses for the distri- 
bution of dense gas in the central kpc of NGC 6946. 



10' 10^ 

Fig. 11.— logioCniiJj'y^CTsus^ogi^oC^hyd/fc) for Go = 100. 
We calculated the error bars by running the la limits of the 
^hcn/^hco+ ratio through the PDR code. 

The dense gas could be concentrated in self-gravitating 
clumps interspersed throughout a less dense molecular 
gas medium. Alternatively, the HCN and HCO+ emis- 
sion could come from a homogeneous component bound 
not by self-gravity, but by the ambient pressure. 

4.1. Self- gravitating Dense Gas and Implications for 
Star Formation 

One hypothesis consistent with our observations is 
that the HCN and HCO+ emission comes from self- 
gravitating dense clumps interspersed throughout less 
dense molecular gas traced by CO. The internal turbu- 
lent pressure of these density peaks is given by: 

Pturb ~ pwfnt (11) 

where w^^t is the internal gas velocity dispersion and p 
is the density of the gas. In these peaks, Pturb is much 
higher than the ambient pressure as measured by Phyd i so 
without self-gravity the density peaks would be unstable. 
This structure is similar in distribution to GMCs in the 
galactic disk, where HCN and HCO"*" are excited in small 
density peaks inside of a envelope traced by CO. We 
have assumed that the turbulent pressure dominates the 
thermal pressure in these peaks, as it does in the star 
forming cores of GMCs. 

We used observations of the nearby Galactic PDR 
in Orion to estimate the Gaussian velocity width 
of the peaks, tJjnt = 2.7 km s"^ (|Boug hton| 1 19781 : 
iTielens fc HollenbachI [l985h . In gSH we found that the 
number density in the inner region of NGC 6946 is con- 
stant with Phyd and i?avg; using nnno ^ 10^ cm^'^ im- 
plies Pturb ^ 10^ cm~'^ K, larger than the calculated 
hydrostatic midplane pressure everywhere in the galaxy. 
Therefore, self-gravity likely plays an important role in 
holding the density peaks together. 

If the density peaks are self-gravitating, we can equate 
the force per unit volume due to turbulent pressure and 
gravity: 

(12) 

where M denotes the mass of the density peak and r is 
the size of the peak. This reduces to the proportionality 
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M oc v^^^r. Using the relation 



niine oc M/r^ oc constant, 



(13) 



we conclude that 



Vint cx r. 



(14) 

Similarly, iRosolowskv &: Blit3 (|2005l ) found for giant 
molecular clouds in the center of M64 that ?;int oc r^-^^^-'^ ^ 
significantly difi^erent from GMCs in the disks of local 
group galaxies. M64 is a galaxy with high ambient pres- 
sure and molecular gas content similar to the bulges of 
normal spirals. The constraint was determined from ob- 
servations of ^^CO, which is better than ^^CO at map- 
ping dense gas due to optical depth effects, but does not 
map the density peaks as well as either HCN or HCO"''. 
The authors then showed that the clouds were overpres- 
surized compared to the hydrostatic midplane pressure 
by calculating the internal pressure from their measure- 
ments and assuming turbulent pressure was dominant. 
Here we find the density peaks are apparently overpres- 
surized and that their number densities are constant, and 
follow the argument in the opposite direction to derive 
the power law they constrained exp e rimen tally. 

In the Milky Way, ISolomon etall (|1987D used CO ob- 
servations of disk GMCs to find that the surface density 
of the cloud is approximately constant with cloud size, 
and therefore the average density is inversely related to 
the cloud size. They constrained the relation between 
the size of the cloud, s, and its velocity dispersion to be 
Vint oc gO-Sio os (note the distinction between size of the 
cloud s and size of the density peak r). Since CO traces 
gas less dense than HCN and HCO"*", these results are 
not necessarily in conflict; density peaks in Milky Way 
clouds may follow the same relation as the density peaks 
in NGC 6946 even though the clouds as a whole do not. 

As i?avg decreases, the rise in both /hcn and /hco+ 
implies that the amount of dense gas available for star 
formation is increasing. Since the density of the peaks is 
constant, this could result from an increase in the size 
of individual density peaks, in the number of density 
peaks above the necessary threshold, or some combina- 
tion thereof. If the effect is due to the number of density 
peaks increasing, this may imp ly a funda mental unit of 
star formation as d iscussed in iGao k, Sol omon (2004bj) 
andl^yu et al.l ()2005f ). However. iKrumhoIz fc Thompson! 
()2007() pointed out a potential hole in this argument: the 
emission of any molecule with a critical density higher 
than the median density is expected to be correlated with 
star formation, regardless of whether the emission and 
star forming regions are directly associated. Further in- 
vestigation would be aided by comparing our integrated 
intensities with the FIR intensity at each point. Unfor- 
tunately, FIR observations with 6" resolution are not yet 
possible, even with the Spitzer space telescope. 

4.2. Homogeneous and Pressure Bound Dense Gas 
Tracers 

Alternatively, the emission from HCN and HCO+ 
could come from homogeneous gas bound by the ambient 
pressure. This situation is similar to what is observed in 
the Milky Way center, where we see emission from HCN 
and CS everywhere we see CO, within the sensitivity lim- 
its (|Helfer fc Blitzlll997al and references therein). That 
is, despite the best efforts to resolve the dense gas tracer 



emission into clumps, it still appears to be diffusely dis- 
tributed. The mean number density of the gas appears 
to be larger than the critical density of the dense gas 
tracers, and there is no need to appeal to self-gravitating 
dense clumps to explain the observed emission. In this 
model, the dense gas filling factor is much larger than 
that of the model described in the previous section. 

To construct a similar situation in the central kpc of 
NGC 6946, support of the molecular gas layer against 
the gravitational potential must come not from turbulent 
pressure, but from thermal pressure: 



Pt 



therm 



= nkT. 



(15) 



The central hydrostatic midplane pressure is 5 x 10^ cm^'^ 
K; if the support comes from turbulence with a charac- 
teristic speed of 3 km s~^, the implied density would be 
~ 2 X 10'* cm"'^, which is too low to excite HCN. In- 
stead, if the thermal pressure provides the support, then 
for kinetic temperatures in the range of 50-200 K, con- 
sistent with those in the center of the Milky Way, the 
derived densities are 2.5 x 10^ — 1 x 10^ cm~^. These 
densities are high enough to excite the dense gas tracers, 
and are marginally consistent with the number densities 
calculated via the FDR code. Temperatures as low as 50 
K are consistent with the FDR calculations with Go ~ 
1000. Higher temperatures would require larger values 
of Go, but the radiation field is not constrained by our 
observations. 

Even at a distance of i?avg ^ 400 pc from the center, 
thermal pressure is likely to play an important role. In 
this region, we know at least some of the gas has n ~ 10^ 
cm~^, because we observe emission from the dense gas 
tracers. Using the FDR code tells us that the tempera- 
ture of this dense gas is ~ 30 K, which implies a thermal 
pressure of 3 x 10^ cm~^ K. The hydrostatic midplane 
pressure for these points is < 5 x 10^ cm~^ K, so thermal 
pressure in this region can provide a significant portion 
of the necessary hydrostatic support, if not all of it. 

Determining whether hydrostatic support at the galac- 
tic center comes from turbulent or thermal pressure 
would point out which of the self-gravitating dense gas 
model or the homogeneous and pressure bound dense gas 
tracer model is correct. We can do this b y constraining 
Ga for the central 10" of NGC 6946 (see iKramerlTall 



2005l) for details on constraining Go); if Go is large 
(~ 10^) the gas is likely warm and thermal pressure will 
dominate, and if it is small (~ 10^) the gas is cool and 
turbulent pressure is more i mportant. Unfortun ately, the 
constraint Go ^ 500 from IStacev et al.l (|1991h was cal- 
culated from C II observations with a resolution of 55" ; 
better resolution than this is necessary to differentiate 
between the two models. 

4.3. Hydrostatic Midplane Pressure and Line 
Intensities 

There is a robust correlation between /hcn and Phyd 
(Fig. [12] {Top)). We fit the relation between the two 
variables with the Monte Carlo routine, which yielded: 



login /i 



10 Jhcn 



-7.0± 1.3+ (1.2±0.2)logi 



hyd 



(16) 

where Phyd is in cm K. It is clear that the total emis- 
sion from HCN increases as Phyd rises. 
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Fig. 12. — Top: Hydrostbyd'niid^?ane pr'esfeure p lott ed against 
^HCN- The solid line is the fit described in eqn. 1161 Bottom: 
Hydrostatic midplane pressure plo tted against 7hco+- The solid 
line is the fit described in eqn. 1171 

We performed the same fit for /hco+ on the coarser 
8.6"x7.5" grid (Fig.MiBottom)). The Monte Carlo fit 
yielded: 



logio ^HCo+ = -4.0 ± 1.0 + (0.7 ± 0.1) log 



10 



Pv 



hyd 



(17) 

where Phyd is in cm ^ K. Again, it is clear that /hco+ 
rises with increasing Phyd- These two relations appear to 
have significantly different slopes, but that may be due 
to the differing resolution of the /hcn and /hco+ maps. 

These fits, along with the radial fits in i i3.3l and the hy- 
drostatic midplane pressure map fPig. llOp . indicate that 
the molecular gas in the central kiloparsec of NGC 6946 
is under high pressure , confirming a prediction made in 
ISpergel fc Blitd (|1992f ). Emission from dense gas trac- 
ers like HCN and HCO+ is widespread, and their inte- 
grated intensity ratios to CO are larger than those fur- 
ther from the center. This gas is in a markedly different 
state than the local interstellar medium in the Galactic 
disk where the pressure is approximately three orders of 
magnitude lower. The structure of the molecular gas dis- 
tribution, whether it is self-gravitating or homogeneously 
distributed and pressure bound, also bears little resem- 
blance to the local interstellar medium. 

4.4. Hydrostatic Midplane Pressure and Molecular Gas 

The hydrostatic midplane pressure has a nearly linear 
relationship with the Sh2/S]hi ratio, Rmoi, over three or- 
ders of magnitude in pressure (jBlitz fc Rosolowskvll200l 



100 




10^ 10' 10° 

p. ./k [cm"^ K] 

Fig. 13. — Hydrostatic rnJifplanc pressure vs. iJmol = /^Hl- 
The HI surface density has been assumed constant over the range 
of this plot. The error bars are determined only from the error in 
Sjjj ; we do not include error in the HI surface density. The line is 
the fitted relation for normal galaxies from BR06. 



|2006| ). Fig. [T2] plots -Rmoi against the hydrostatic mid- 
plane pressure for the points in the inner central 350 pc 
of NGC 6946, along with the fit to normal galaxies from 
BR06. The fit follows the points reasonably well, though 
the data are slightly displaced toward higher Phyd- How- 
ever, BR06 did see comparable scatter in their original 
sample around this relation. This paper supplements the 
previous work by adding more points at the high pres- 
sure (10^-10^ cm~'^ K) end of the relationship, albeit 
with the assumption that the HI surface density is con- 
stant. With such large molecular gas surface densities, 
the HI gas should not contribute a significant amount to 
the total gas surface density anyway. 

Both Phyd and Pmoi depend on /co, and plotting one 
as a function of the other can lead to misleading cor- 
relations. However, the fit found in BR06 was deter- 
mined from many points with Pmoi ~ 1, and the authors 
showed that the relationship was physically significant in 
spite of the self-correlation issues. We include this plot 
to demonstrate the continued applicability of the fit in 
the high pressure regime. 

5. CONCLUSIONS 

Using the BIMA and CARMA interferometers, we pro- 
duced high-resolution maps of the dense gas tracers HCN 
and HCO+ in NGC 6946. The maps are rich with detail, 
and we have 13 independent data points across the inner 
^750 pc of the galactic center in HCN and CO, and 7 
in HCO+. The integrated intensities of HCN and HCO+ 
peak near the center of the galaxy, and fall off from there. 

Using BIMA SONG CO 1-0 and 2MASS Ks band im- 
ages, we calculated an approximation to the hydrostatic 
midplane pressure in the center of NGC 6946. We used a 
PDR model to calculate the number density in the dense 
gas regions traced by HCN and HCO+, and found that 
the number density in these regions does not vary as a 
function of radius or the surrounding hydrostatic mid- 
plane pressure. We explored two hypotheses for the dis- 
tribution of the dense gas in the central kpc of NGC 6946. 
If the dense gas is concentrated in clumps, we showed 
that self-gravity implies that Wjnt cx: r. However, we also 
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explored another plausible dense gas distribution; the 
dense gas could instead be homogeneously distributed 
throughout the galactic center if thermal pressure pro- 
vides enough support for hydrostatic equilibrium. We 
then demonstrated clear correlations between the HCN 
and HCO"*" integrated intensities and the midplane pres- 
sure. Finally, we confirmed the validity of the relation 
between hydrostatic midplane pressure and the molecu- 
lar to atomic gas ratio found in BROS in the high pressure 
regime. 
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